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We study the production of Mueller-Navelet jets at hadron colliders in the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) framework. We show that a measurement of the relative azimuthal angle A$ 
between the jets can provide a good testing ground for corrections due to next-leading logarithms 
(NLL) . Besides the well-known azimuthal decorrelation with increasing rapidity interval An between 
the jets, we propose to also measure this effect as a function of R = fe/fei, the ratio between the 
jets transverse momenta. Using renormalisation-group improved NLL kernel, we obtain predictions 
for da I dAndRdAQ . We analyse NLL-scheme and renormalisation-scale uncertainties, and energy- 
momentum conservation effects, in order to motivate a measurement at the Tevatron and the LHC. 
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O . I. INTRODUCTION 

Mueller-Navelet jet production [l[ in hadron-hadron scattering is a process in which a jet is detected in each of 
the forward directions with respect to the incident hadrons. This process is characterized by two hard scales: k\ 
and &2, the transverse momenta of the forward jets. When the total energy of the collision y/s is sufficiently large, 
corresponding to a large rapidity interval between the jets Ary^lnQi/fciA^), Mueller-Navelet jet production is relevant 
for testing the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach [2]. 

In fixed-order perturbative QCD calculations, the hard cross section is computed at fixed order with respect to a s . 
The large logarithms coming from the strong ordering between the hadrons scale and the jets transverse momenta 
are resummed using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation Q for the parton 
i' densities. However in the high-energy regime, other large logarithms arise in the hard cross section itself, due to the 
strong ordering between the energy yfs and the hard scales. These can be resummed using the BFKL equation, at 
i " leading (LL) and next-leading (NLL) logarithmic accuracy [1, H[ ■ 

On the phcnomcnological side, a first attempt to look for BFKL effects was performed at the Tevatron (Run 1), 
using measurements of cross-section ratios (for same jet kinematics and two different center-of-mass energies squared 
s and s) that are independent of the parton densities and allow to study more quantitatively the influence of the 
high-energy effects. The data [f| overestimate the LL-BFKL prediction (s / s) 4q ln ( 2 ), however it has been argued @ 
~-r ' that the measurement was biased by the use of upper Et— cuts, the choice of equal lower Et— cuts, and hadronization 
, corrections. As a result, these tests on the relevance of the BFKL dynamics were not conclusive. 

On the theoretical side, it was known that NLL corrections to the LL-BFKL predictions could be large due to the 
appearance of spurious singularities in contradiction with renormalization-group requirements. However it has been 
realised 0, [§[ that a renormalisation-group improved NLL regularisation can solve the singularity problem and lead 
to reasonable NLL-BFKL kernels (see also @ for different approaches) . This motivates the present phenomenological 
• • ■ study of NLL-BFKL effects in Mueller-Navelet jet production. Our analysis allows to study the NLL-BFKL framework, 
and the ambiguity corresponding to the dependence on the specific regularisation scheme. Our goal is to motivate 
. further measurements at the Tevatron (Run 2) and at the LHC. 

In Ref. [llj and [ll|, such phenomenological investigations have been devoted to the proton structure function 
and forward-jet production in deep inelastic scattering. The NLL-BFKL effects were taken into account through an 
"effective kernel" (introduced in Q) using three different schemes (denoted S3 and S4 from Q and CCS from 0]). 
While for the structure function analysis the NLL corrections didn't really improve the BFKL description, it was 
definitively the case in the forward-jet analysis. 

The present study is devoted to the A<I> spectrum, where A$ is the relative azimuthal angle between the Mueller- 
Navelet jets. We implement the NLL-BFKL effects following [l(| and [llj], using the S3 and S4 schemes. We study 
the magnitude of the NLL corrections with respect to the LL-BFKL results. We confirm the expectations [ll| that 
those corrections slow down the azimuthal decorrelation with increasing Arj. 
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We propose to also investigate this effect as a function of R = k%/ki, the ratio between the jets transverse momenta. 
This is inspired by the results of [Til ] which showed that NLL-BFKL corrections have more impact on the forward-jet 
cross-section when the measurement is sensitive to different values of (the forward-jet equivalent of) R. We obtain 
predictions for da h J J / dAi~idRdA$> and show that this would allow for a detailed study of the NLL-BFKL approach 
and the QCD dynamics of Mueller-Navelet jets. 

The plan of the paper is the following. In section II, we present the phenomenological NLL-BFKL formulation of 
the Mueller-Navelet jet cross-section. In section III, we introduce the observable da h J / dA-qdRdAQ relevant to 
study the A$ spectrum. In section IV, we present the predictions obtained using the S3 and S4 schemes and compare 
them with LL-BFKL predictions. We also discuss the dependence of our results with respect to the choice of the 
renormalization scale determining a s , and we estimate the impact of energy- momentum conservation effects. Section 
V is devoted to conclusions and outlook. 



II. MUELLER-NAVELET JETS IN THE NLL-BFKL FRAMEWORK 



Mueller-Navelet jet production in a hadron-hadron collision is represented in Fig.l with the different kinematic 
variables. We denote ^/s the total energy of the collision, k\ and fc 2 the transverse momenta of the two forward jets 
and x\ and cc 2 their longitudinal fraction of momentum with respect to the incident hadrons as indicated on the 
figure. A<I> = 7r— 4>i+4>2 measures the relative azimuthal angle between the two jets, as <pi and 4>2 are the jets angles 
in the plane transerve to the collision axis. In the following, we consider the high-energy regime in which the rapidity 
interval between the two jets Ar) = log(ccia; 2 s/fci/c 2 ) is assumed to be very large. Following the phenomenological 
NLL-BFKL analysis of [13, ED] , one obtains the Mueller-Navelet jet cross section: 



da hf ^ JXJ _ a.(fc?)a.(*2) , ^ ^ f ^ ^ f d 1 (kf 

U2 J)„2, 



f eff (x 1 ,kl)f eff (x 2l kl) jr e^ k ^ff^^ k ^ + ^ (1) 



dxidx 2 dkfdk 2 dA^ ikfk^ 
with the complex integral running along the imaginary axis from 1/2 — zoo to 1/2+zoo. The running coupling is 

a(k*)=a s (e)N c /n=[b\og(e/Al CD )r 1 , b= ■ (2) 

Let us give some more details on formula ([T]). 

• The NLL-BFKL effects are phenomenologically taken into account by the effective kernels XeffiPiJ-,®)- For 
p = 0, the scheme-dependent NLL-BFKL kernels provided by the regularisation procedure xnll{i,u) depend 
on 7, the Mellin variable conjugate to k\jk\ and u>, the Mellin variable conjugate to s/sq where sq — k\k2 is 
the energy scale. In each case, the NLL kernels obey a consistency condition 7] which allows to reformulate the 
problem in terms of Xeff{l,OL). The effective kernel Xeffil, &) is obtained from the NLL kernel xnll{i,w) by 
solving the implicit equation Xeff — Xnll{i, ol Xeff) as a solution of the consistency condition. 

In the case of the S3 and S4 schemes 0] (in which xnll is supplemented by an explicit a dependence), we will 
extend the regularisation procedure to non zero conformal spins and obtain Xnll(p, 7> <^); this is done in the 
Appendix. Then the effective kernels XeffijPt 7j ot) are obtained from the NLL kernel by solving the implicit 
equation: 

Xeff = Xnll(p, 7, a Xeff) ■ (3) 

• In formula fTJ), the renormalisation scale determinig a is k 2 = k\k2, in agreement with the energy scale sq [I3.[T3|. 
In Section IV, we shall test the sensitivity of our results when using k 2 — A fcife and varying A. This is done 
using formula JT]) with the appropriate substitution [llj 

a(fci/c 2 )->a(Afcifc 2 ) + & a 2 (fci/c 2 )log(A) , (4) 

and with the effective kernel modified accordingly following formula ([3|). We also modify the energy scale into 
,s = A kik 2 . 

• It is important to note that in formula ([T]), we used the leading-order (Mellin-transformed) impact factors. We 
point out that the next-leading impact factors are known [151] . and that in principle, a full NLL analysis of 
Mueller-Navelet jets is feasible, but this goes beyond the scope of our study. Also, our formula is different from 
the one proposed in [l6j], because the authors considered the cross-section integrated with respect to the jets 
transverse momenta. This leads to a modification of the jet impact factors which results in an extra factor 
7 _1 (1 — 7) _1 in the integrand of (p}. Also it modifies the effective kernel (see 16] where the S3 scheme was 
considered). 



fa, yi= hx(xiy/s/ki) 




Arj = h-i(x 1 x 2 s/(fak 2 )) 



fa, 2/2 = - ln{x 2 ^/s / k 2 ) 



FIG. 1: Mueller-Navelet jet production in a hadron-hadron collision. The kinematic variables of the problem are displayed, s 
is the total energy squared, k\ (yi) and 2v2 (1/2) are the transverse momenta (rapidities) of the jets and x\ and xi are their 
longitudinal momentum fraction with respect to the incident hadrons. At] is the rapidity interval between the hard probes. 



In formula ([T]), f e ff(x,k 2 ) is the effective parton distribution function and resums the leading logarithms 
log(k 2 / Aq CD ). It has the following expression 



f eff (x, k 2 ) = g(x, k 2 ) + ^ (q(x, k 2 ) + q(x, k 2 )) 



(5) 



where g (resp. q, q) is the gluon (resp. quark, antiquark) distribution function in the incident proton. Since the 
Mueller-Navelet jet measurement involves perturbative values of k\ and k 2 and moderate values of X\ and x 2 , 
formula ^ features the collinear factorization of f e ff, with kf and k 2 chosen as factorization scales. 



By comparison, the LL-BFKL formula is formally the same as ((T]), with the substitutions 



Xeffip, 7, a) -> Xll{p, 7) = 2V(1) - V fl - 7 + y) ~ * ( 7 + ^f) ' " (fc2 
where ^(7) =d log T ("yj/dy is the logarithmic derivative of the Gamma function. 



a — const, parameter 



(6) 



III. THE A$ SPECTRUM 



We would like to study the azimuthal decorrelation of the Mueller-Navelet jets as a function of their transverse 
momenta k\ and k 2 and rapidities y\ and y 2 : 



2/i = log 



2/2 = - log 



X 2 y/S 



k\ J V fa 

Let us first introduce kinematic variables suitable for our problem: we change the variables in ([T]) to the variables 



(7) 



y = 



2/1 + 2/2 



Q = y/fafa 



and R = 



fa 



(8) 



One obtains 



da 



hh^JXJ 



dATjdydQdRdA® 



a s (Q 2 /R)a s (Q 2 R) 2 2 
Q3 xif eff {x u Q /R)x 2 f e ff{x2,Q R) 

00 /-1/2+0O 7 

/ ZI^- 2 7 e a(Q 2 )xe//b,7,a(Q 2 )]A,,+ipA* 



(9) 
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FIG. 2: Left plots: values of d- p (Aij,R = 1) (see formula (|13[0 entering into the A$ spectrum for the rapidity intervals 
Ay — 6, 8, 10; upper plot: LL-BFKL, middle plot: S3 scheme, lower plot: S4 scheme. Right plots: ratios dp /ap] upper 
plot: S3 scheme/LL-BFKL, lower plot: S4 scheme/LL-BFKL. 



We are interested in the following observable, suitable to study the azimuthal dccorrelation of the jets as a function 
of their rapidity separation Arj and of the ratio of their transverse momenta R : 



2- 



d,a 



da 



= 1 



cr p (Af7,i?)cos(>A$) . 



dA-qdRdA^ I dA-qdR 
We have expressed the normalized cross-section (1101) in terms of the Fourier coefficients 

da \ 1 f , . T , . T , da a p (Ar],R) 



(cos(pA$)) = 



dArjdR 



dA$ cos(pA$) 



dArjdRdA® a {Arj,R) 



(10) 



(11) 



with the cross-sections <7 P (A?7, R) obtained from ^ and given by 



a p (Ar),R) 



dQ 

Q 3 



a s (Q 2 /R)a s (Q 2 R) 



dy Xif e ff{x 1 ,Q 2 /R)x 2 feff(x2,Q 2 R) 
/■l/2+oo j 

/ -Lr- 2 ^ e a(Q 2 )Xe//b,7,a(Q 2 )]Ar, 

il/2-oo 2i7T 



(12) 



The kincmatical cuts Q > Et and y < < y < y > for the Q and y integrations in (112jl will be specified later, when we 
discuss the Tevatron and LHC kincmatical ranges. 

For the sake of comparison between BFKL LL and NLL effects, we define the following quantities, free of parton 
distribution functions: 



a p (Ar), R, a) 



l/2+oo , 
l/2-oa 2l7T 



(13) 



Note that in the LL-BFKL case in which a does not depend on Q 2 , one has a p /ao — a v ja^. We shall compare the 
LL and NLL values of <7 P (A?7, R, 0.16) for R = 1 and A-q = 6, 8, 10. The comparison is shown on Fig. 2 where we 
consider both the S3 and S4 NLL schemes. 

The cross sections a p are displayed as a function of p and, as expected for the rather large values of Arj considered, 
we see that do is the largest cross section, and its increase with rapidity is stronger at LL compared to NLL. For 
pj^O, a p decreases as a function of Arj, and the ratios cr p LL ja p L between the NLL and LL contributions show that 
the decrease is faster at NLL except for p = 1 and p ~ 2 (and for p = 3 the rapidity dependences at LL and NLL are 
comparable). 




FIG. 3: The Mueller-Navelet jet A"l> distribution (|10|l for Tevatron (run 2) kinematics in the BFKL framework at LL (upper 
plots) and NLL-S4 (lower plots) accuracy. Left plots: R = 1 and Ar/ = 6, 7, 8. Right plots: Ar/ = 8 and R — 1, 1.5, 2. 



IV. RESULTS FOR MUELLER-NAVELET JET A$ DISTRIBUTIONS 

In this section, we show the results for the A$ distribution obtained with formulae (fT0|) and (TT2|) . As shown in 
Fig. 2, dp decreases as a function of p, and the decrease is faster at NLL compared to LL (and is similar for both 
schemes S3 and S4). As a result, including 20 terms in the sum over p in (fit))) is enough in the S3 and S4 cases. 
However at LL, one has to include more terms depending on the value of A77 and R. 

We choose to apply the rapidity cut \y\ < 0.5 which enforces a symmetric situation y^ ~ — y\. For the transverse 
momentum cut Et, we will consider two options corresponding to the Tevatron and the LHC possibilities in terms 
of kinematical reach: Et = 20 GeV for the Tevatron (Run 2) and Et = 50 GeV for the LHC. We recall that the 
respective center-of-mass energies are a/s = 1960 GeV and a/s=14 TeV. 

We point out that our NLL-BFKL predictions for the observable (jTU)) are parameter free. In the LL-BFKL case 
that we consider for comparisons, the only parameter a is fixed to the value 0.16 obtained in [l7j by fitting on 
forward-jet data from HERA. By contrast, in the NLL-BFKL case, the value of a is imposed by the renormalisation 
group equations. 

A. Comparison between LL and NLL BFKL predictions at the Tevatron and the LHC 

In Fig. 3 and Fig. 4, we display the observable (fT0|) as a function of A$, for Tevatron and LHC kinematics respectively. 
The results are displayed for different values of A77 and R and at both LL and NLL accuracy (in this case, the S4 
scheme is used). In general, the A$ spectra are peaked around A$ = 0, which is indicative of jet emissions occuring 
back-to-back. In addition the A$ distribution flattens with increasing Ar] = yi — y2 or with R=k2/ki deviating from 
1. Note the change of scale on the vertical axis which indicates the magnitude of the NLL corrections with respect 
to the LL-BFKL results. The NLL corrections slow down the azimuthal angle decorrelations for both increasing A77 
and R deviating from 1. 

In the BFKL framework, the A<I> dependence of the spectrum (fT0|) is larger at NLL than at LL. However, this 
A$ dependence is still smaller than in the fixed-order pQCD approach, in which the back-to-back peak is quite 
pronounced. Therefore a measurement of the cross-section da hh ^ JXJ /dAr/dRdAQ at the Tevatron (Run 2) or the 
LHC would allow for a detailed study of the QCD dynamics of Mueller-Navelet jets. In particular, measurements 
with values of Ar] reaching 8 or 10 will be of great interest, as these could allow to distinguish between BFKL and 
DGLAP resummation effects and would provide important tests for the relevance of the BFKL formalism. In addition, 
measuring the normalized cross-section (|10[) could help reducing the biases which altered previous measurements 0, @] . 




FIG. 4: The Mueller-Navelet jet A"I> distribution (fT0)l for LHC kinematics in the BFKL framework at LL (upper plots) and 
NLL-S4 (lower plots) accuracy. Left plots: R = l and A 77 = 6, 8, 10. Right plots: Ar? = 10 and R = 1, 1.5, 2. 



The DO collaboration at the Tevatron (Run 1) did measure the azimuthal angle distribution between two jets 
[l8|, but they were not separated in rapidity by more than 5 units, in which case we do not expect the BFKL 
predictions to be relevant. Nevertheless, fixed order QCD predictions at next-to-leading order failed to describe the 
data, underestimating the decorrelation. In contrast, NLL-BFKL calculations overestimate the decorrelation [ll|. 
Solving this puzzle likely requires to measure Mueller-Navelet jets with higher values of A77. 

B. Scheme and scale dependence 

Our previous results in the NLL-BFKL case were obtained with the S4 scheme. As shown in Fig. 2, the S3 scheme 
leads to similar results for the quantities a p (Ar),R) and this is also true for the cross-sections a p (Ar),R) (formula 
(fT2"|) ) that actually enter in the formulation of the observable (fTU|) . There are some differences between the S3 and 
S4 scheme, but they tend to cancel when computing the ratios a p /ao to obtain the A<I> spectrum. Therefore the 
results obtained with both schemes are almost not distinguishible, as displayed on the left plots of Fig. 5. Let us also 
point out that the pdf uncertainties cancel in the same way, and that the effects (not implemented here) due to the 
next-to-leading order jet impact factors would be suppressed too. 

Let us now study the renormalization scale dependence of the NLL-BFKL description of Mueller-Navelet jets. 
Previously, the choice was k\ki = Q 2 and we now test the sensitivity of our results when using Q 2 /2, and 2Q 2 . We use 
formula |T]) with the appropriate substitution a(Q 2 ) — > a(XQ 2 )+b a 2 (Q 2 ) log(A) and with the effective kernel modified 
accordingly following formula ^ . We also modify the energy scale Q 2 — > A Q 2 . The results are shown on the right 
plots of Fig. 5, and the dependence on the choice of scale turns out to be quite small, of about 5 percent, except for 
A$ close to 0, in which case the uncertainty reaches 20 percent. 

C. Energy- momentum conservation effects 

The analytic expression of the BFKL cross-section ([T]) lacks energy-momentum conservation, because these effects 
are formally higher-order corrections in this framework. However it has been argued pjj [2(| that of the terms 
which conserve energy- momentum could be numerically important for phcnomcnological analysis. Therefore we shall 
estimate their magnitude for the observable (fT0|) . In order to do so, we will use the proposal of which amounts 
to substitute Arj in (fTJ) by an effective rapidity interval y e ff- More advanced Monte Carlo approaches were later 
developed [20j, but we choose to stick to more insightful analytic calculations. 




FIG. 5: Resumation-scheme and renormalization-scale dependencies of the Mueller-Navelet jet A$ distribution (|10p in the 
NLL-BFKL framework. Upper plots: R = l, A?7 = 8 and Tevatron (run 2) kinematics; lower plots: R—l, Ar/ — 10 and LHC 
kinematics. The left plots show a comparison of the S3 and S4 schemes while the right plots display results obtained with the 
three renormalization scales Q 2 /2, Q 2 , 2Q 2 . 




FIG. 6: Effects of energy conservation on the Mueller-Navelet jet A$ distribution for Ar/ — 10 and LHC kinematics. Left plot: 
R—l; the effect is minimal. Central plot: R — 1.1, right plot: R — 1.3; the azimuthal correlation increases with R deviating 
from 1 (instead of decresing) after energy-momentum conservation is included. 



The effective rapidity is defined in the following way 

/ r da°^ \ 1 f d(T LL - BFKL 

y eff (p,Q,R,Arj,y) = Arj \J # coefefl dArjdydQdRdA $ j J d ^ os ^ dA7]dydQdRdA ^ ( 14 ) 

where da ^ 01 ^ is the exact 2 — > 3 contribution to the Wi— > JXJ cross-section at order [2l[, and da LL ~ BFKL is 
the LL-BFKL result. One has y e ff(Ar]— > oo) = Ar]. In this way, when used in the expansion of the cross-section 

with respect to a s is exact up to order a, while the large Arj limit is unchanged. To compute da°( a *\ we used the 
standard jet cone size R cu t = 0.5 when integrating over the third particle's momentum. The main feature of y e ff is 
that it is only slightly smaller than Arj for R = 1 , but that it decreases quickly with R deviating from 1 [l9[ . 
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FIG. 7: The Mueller-Navelet jet A$ distribution JlO} for CDF kinematics and R = 1. Left plot: NLL-BFKL predictions for 
A?7 = 6, 8, 10, 11. Right plot: comparison with the LL-BFKL result and calculations taking into account energy conservation, 
this effect is small as R = 1. 



As shown in Fig. 6, where the observable (fTT)|) is plotted for LHC kinematics and Ar/= 10, this behavior is confirmed. 
Indeed, when R = 1 the effect is minimal, the azimuthal correlation is only slightly bigger with energy momentum 
conservation. By contrast when R^l, the azimuthal correlation is much bigger with energy momentum conservation 
than without, and the effect is more and more important as R deviates from 1. Therefore the modification of the 
Acf> spectrum with respect to R is a measure of the role of energy-momentum conservation effects: without them the 
azimuthal correlation decreases with R deviating from 1 while it is the opposite if such effects are included. 



D. Mueller-Navelet jets at CDF 



The CDF collaboration recently installed detectors called Miniplugs in the forward and backward regions. These 
detectors allow to increase the acceptance in rapidity and transverse momentum to measure very forward jets. It will 
be possible to measure jets separated in rapidity by more than 10 units and with transverse momenta as low as 5 GeV. 
It is also worth pointing out that while the CDF Miniplug detectors are not prefectly suited for energy measurements 
(the jet containment will be poor: the depth of the calorimeters is only one A), they are especially interesting in the 
case of the observable studied here, which focuses on the difference in azimuthal angle between the jets. 

The NLL-BFKL predictions for the Mueller-Navelet jet A<& distribution with CDF kinematics is represented in 
Fig. 7. With such low values of transverse momenta and large values of rapidity interval between the two jets, it is 
also likely that saturation effects could play an important role. First estimations [22j (obtained with less favorable 
kinematics) indicate so when considering saturation effets damping the LL-BFKL exponential growth. Stud yin g 
saturation effects with NLL-BFKL growth certainly deserves more study. First steps have been taken in Ref. [23j |. 
but the problem of phenomenology for hadron colliders has yet to be addressed. 



V. CONCLUSION AND OUTLOOK 



We have investigated the decorrelation of Mueller-Navelet jets with respect to their relative azimuthal angle A$ in 
the BFKL framework at NLL accuracy. Using renormalization- group improved NLL kernels xnll{p, 7, in the S3 
and S4 schemes, the NLL-BFKL effects were taken into account through an effective kernel obtained from the implicit 
equation ([3]) . This allowed our phenomenological study of NLL-BFKL effects in Mueller-Navelet jet production. Our 
present goal is to motivate future measurements at the Tevatron (Run 2) and at the LHC |25j. A future comparison 
with the data will require to adapt our predictions to experimental cuts and perhaps to less differential cross sections. 



The present study, devoted to the A<I> spectrum (|10p . confirms the expectations that when increasing the rapidity 
interval between the jets A?7, the decorrelation increases, and that NLL corrections decrease the azimuthal decorre- 
lation with respect to the LL-BFKL results. We also investigated this effect as a function of R = k^jkx, the ratio 
between the jets transverse momenta: when R deviates from 1, the azimuthal decorrelation increases. Our predictions 
were obtained with standard expectations of Tevatron and LHC kinematical possibilities. However, we also presented 
predictions for the Mueller-Navelet jet A<i> distribution having in mind the CDF forward detector which features a 
quite favorable kinematical reach (Q>5 GeV and A?7>10). 

For the observable (fit))) , we noticed that the differences between the different schemes are quite small, while the 
dependence on the choice of renormalization scale is of about 5 percent in general and reaches 20 percent around 
A$ = 0. Energy- momentum conservation effects are minimal for R=l, but they increase quite rapidly as R deviates 
from 1. In fact, they reverse the trend discussed above: with energy-momentum conservation implemented, the 
azimuthal decorrelation decreases as R deviates from 1. 

With such low values of transverse momenta and large values of rapidity interval, Mueller-Navelet jet measurements 
would allow for a detailed study of the QCD dynamics of Mueller-Navelet jets, both for investigating fixed-order pQCD 
versus BFKL predictions, but also with respect to possible saturation effects. In these contexts, the measurement 
of the A<i> integrated cross-section would be very interesting by itself, but a realistic phenomenological study should 
incorporate the next-to-leading order jet impact factors in the calculation. Indeed, their effect will not be suppressed 
as it likely is in the case of the normalized cross-section we have studied in this paper. 



Acknowledgments 

We would like to thank Robi Peschanski for commenting on the manuscript. CM. is supported in part by RIKEN, 
Brookhaven National Laboratory and the U.S. Department of Energy [DE-AC02-98CH10886]. 



Appendix: The S3 and S4 schemes for non-zero conformal spins 



In this Appendix, we show how to extend the regularisation procedure of [Tj to non zero conformal spins p ^ 0. We 
obtain xnll(p,J,u) for the S3 and S4 schemes (recently two preprints appeared where the S3 scheme [16[ and the 
other Salam schemes [H[ have also been extended). 

The starting point is the scale invariant (and 7<->l — 7 symmetic) part of the NLL-BFKL kernel 



Xi (P, 7)= 2^(3) + 



1 + 56 C(2) 



Xll{p,i)~ ^xl L (P,l) + l [^"(7 + f) + <(l-7 + f) 



with b given in 



"2 [0(P>7) +<Hp> !-7)] 
Xll given in ((6]), and 



7T 2 C0S(7T7) 



4 sin 2 (^7) (1 - 27) 



N f \ 2 + 37(1-7) 
WJ (3-27)(l + 2 7 ) 



- 1 



7(1 - 7) 



4>{p, 7) = 



(-1) 



k=0 

k r 



k + j + p/2 



p- 



ip'(k + l)-ip'(k+p 
,,fk + p + l 



-1) 



N f 



A c V2(3-2 7 )(l + 2 7 ) 



ip(k + p + 1) - j>(k + I) 
fc + 7+p/2 



#0p 
$2p 



(15) 



(16) 



Note that for the terms on the first line of (fTI)]) inside the curly brakets, we have corrected the signs with respect 
to Ref. [Ill, where they are misprinted (the signs are correct in Ref. [271]). As is the case for xll(p,j), the kernel 
Xi(p, 7) has poles at 7 = — p/2 and 7=1 +p/2. The pole structure at 7 = — p/2 (and by symmetry at 7 = 1 +p/2) 
is: 



XiiP, 7) 



1 



diip) 



2(7 +§) 3 (7+i) 2 



•0(1) 



(17) 



with 



di(p) 



^-+6[l/>(P + l)-lKl)] + s 



-if/ 



■4^'(p + l) 



iV? / 36 V iV? / 1800 



(18) 



and 



d 2 (p) 



2 2 



[^(p + l)-V(l)]- H + 2 



f \ °o P 



N? y 12 



AT? y 60 



(19) 



Note that Xi(2, 7) also has a pole at 7 = with residue (1 + Nf/N^)/2A. This manifestation of the non-analyticity 
[26| of Xi(.Pj7) with respect to the conformal spin does not alter the stability of the NLL prediction and a careful 
treatment of this singularity is not required. 



A. Extension of the S3 scheme 



The S3-scheme kernel XS3(P, 7,w) is given by 



XS3{p,J,u) = [1 - aA{p)} 



20(1) -0 7 + 



p + 2aB(p) +uj 



-0 1-7 + 



p + 2aB(p) + oj 



+« \ Xi(p,l) + A(p)xll(p,i) + [B(p) + 



Xll{p,i) 



^( 7 + |) + ^( 1 - 7 + |) 



(20) 



with A(p) and B{p) chosen to cancel the singularities of Xi(p, 7) at 7 = — p/2 



A(p) = -di(p)-^(p + l) , S(p)=-da(p) + 2^(p + l)-^(l) 



(21) 



B. Extention of the S4 scheme 



The S4-scheme kernel xsi(p, 7j^) is given by 



Xsi{p,l,u) = xll{p,i) ~ /(P,7) + [1- q; ^(p)]/(p + w + 2aB(p),j) 



^ Xi (p, 7) + A(p)/(p, 7) + B(p) + 



Xll{p,i) 



7+ 2 



1-7 + f 



(22) 



with 



/(P.7) = -^ + — 1 



7+| • 



In this scheme, and i?(p) are given by: 



(23) 



A(p) = -d.ip) - \ 



+ 1) - V'(i) + 



(p+1) 2 



B(p) = -d 2 (p) + -[0(p + l)-V(l)] 



(24) 
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